Novel Attractive Force Between Ions in Quantum Plasmas 
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We report new attractive force between ions that are shielded by degenerate electrons in quantum 
plasmas. Specifically, we show that the electric potential around an isolated ion has a hard core 
negative part that resembles the Lennard- Jones (LJ)-type potential. Physically, the new electric 
potential is attributed to the consideration of the quantum statistical pressure, the quantum Bohm 
potential, as well as the electron exchange and electron correlations due to electron- 1/2 spin within 
the framework of the quantum hydrodynamical description of quantum plasmas. The shape of the 
attractive potential is determined by the ratio between the Bohr radius and the Wigner-Seitz radius 
of degenerate electrons. The existence of the hard core negative potential will be responsible for 
the attraction of ions forming lattices and atoms/molecules, as well as for critical points and phase 
transitions in quantum plasmas at nanoscales. 

PACS numbers: 71.10.Ca, 63.10.+a, 67.10.Hk 



A study of the potential distribution around a test 
charged object is of fundamental importance in many 
physical systems (e.g. dilute charge stabilized colloidal 
suspensions such as latex spheres in water, micelles and 
microemulsions, condensed matters with strongly corre- 
lated electrons and holes, strongly coupled dusty and 
quantum plasmas in laboratory and astrophysical envi- 
ronments, etc.), since its knowledge predicts how a cloud 
of opposite polarity charges will shield a test charge par- 
ticle over a certain radius, which is known as the Debye- 
Hiickel radius in the context of electrolytes and plas- 
mas and the Thomas-Fermi/ Yukawa radius in the con- 
text of condensed matters. The traditional repulsive 
screened Coulomb potential assumes the form cf>(r) — 
(Q/r) exp(— r/A), where Q is the test charge and A the 
shielding radius of the sphere. In a classical electron- 
ion plasma 0, HJ, an isolated ion is shielded by non- 
degenerate Boltzmann distributed electrons, and hence 
one [3j replaces Q by Zie and A by the electron Debye ra- 
dius A_De = (fcsT e /47m e 2 ) 1 / 2 , where Zi is the ion charge 
state, e the magnitude of the electron charge, ks the 
Boltzmann constant, T e the electron temperature, and hq 
the unperturbed electron number density. For a slowly 
moving test charge in collisionless Q and collisional 0, H| 
plasmas, there appear additional far-field potentials de- 
creasing as inverse cube and inverse square of the dis- 
tance between the test charge and the observer. In a col- 
lisionless dusty plasma [749} with Boltzmann distributed 
electrons and ions, a micron-sized negatively charged iso- 
lated dust is shielded by both positive ions and electrons. 
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Here Q equals Z^e and A is replaced by the effective dusty 
plasma Debye radius Xd — A_d c Az>;/ V A|> e + where 
Zd is the number of electrons residing on a dust grain, 
Arji = (ksTi/ iirrnoe 2 ) 1 / 2 the ion Debye radius, Tj the 
ion temperature, n^o — tiq + ZdHdo the ion number den- 
sity, and rido the dust number density. In dusty plas- 
mas with T e ^> Ti, we have Arj rs Arj,;. Furthermore, in 
dense Thomas-Fermi plasmas with an impurity ion (with 
the charge state -Z*) is shielded by non-relativistic de- 
generate electrons, so that Q = Z*e and A is replaced 
by the Thomas-Fermi radius Af = (£f /^Trnoe 2 ) 1 / 2 , 
where the electron Thomas-Fermi energy is denoted by 
£p = (?i 2 /2m*fcB)(37r 2 ) 2 / 3 nQ^ 3 , h is the Planck constant 
divided by 2n, and m» the effective mass of the electrons 
(for example, for semiconductor quantum wells, we typ- 
ically have to* = 0.067m e , where m e is the rest mass of 
the electrons). 

In the past, tremendous progress has been made in 
carrying out systematic theoretical and numerical stud- 
ies of phase diagrams [TiTfr7l | for colloidal systems, dusty 
plasmas, and strongly interacting matter by supposing 
that like-charged particles repel each other due to the 
Debye-Hiickel/ Yukawa repulsive force. However, besides 
the repulsive interaction, there are also attractive forces 
0, between two like charged particles due to the over- 
lapping Debye spheres [l8| and due to the polarization of 
charged particulates by the sheath electric field [19j, l2fj . 
Henceforth, both short range repulsive and long-range at- 
tractive potentials play a decisive role in the theory and 
experiments of phase transitions in colloidal and dusty 
plasmas. 

In this Letter, we present a new attractive force be- 
tween ions that are shielded by the degenerate electrons 
in strongly coupled quantum plasmas that are ubiquitous 
in a variety of physical environments (e.g. the cores of 
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Jupiter and white dwarf stars [2lj, |22j , warm dense mat- 
ters [23]]), in compressed plasmas produced by intense 
laser beams (24|, as well as in the processing devices for 
modern high-technology (e.g. semiconductors [25], thin 
films and nano-metallic structures [26j], etc.) Specifically, 
we shall use here the generalized quantum hydrodynam- 
ical (G-QHD) equations [2(| for non-relativistic degen- 
erate electron fluids supplemented by Poisson's equa- 
tion. The G-QHD model includes the quantum statis- 
tical pressure and the Bohm potential [27H29j |. as well 
as the electron exchange and electron correlation effects. 
We demonstrate here that the electric potential around 
an isolated ion in quantum plasmas has a new distribu- 
tion, the profile of which in special cases resembles the 
Lennard- Jones-type potential. It emerges that the newly 
found electric potential, arising from the static electron 
dielectric constant that includes the combined effects of 
the density perturbations associated with the quantum 
statistical pressure and the quantum force [30h32| ] involv- 
ing the overlapping of the electron wave function due to 
the Heisenberg uncertainty and Pauli's exclusion princi- 
ples, as well as the electron exchange and electron cor- 
relation effects 33] due to electron-I/2 spin, embodies 
a short-range negative hard core electric potential. The 
latter will be responsible for Coulomb ion crystallization 
and oscillations of ion lattices under the new force associ- 
ated with our exponential oscillating-screened Coulomb 
potential in strongly coupled quantum plasmas. 

Let us consider a quantum plasma in the presence 
of non-relativistic degenerate electron fluids and mildly 
coupled ions that are immobile and form the neutral- 
izing background. In our quantum plasma, the elec- 
tron and ion coupling parameters are T e = e 2 / a e k B Tp 
and Ti — Z 2 e 2 /a,ik B Ti, respectively, where a e ~ a% — 
(3/471-no) 1 / 3 is the average interparticle distance, and 
Tp = {Ti 2 /2m if k B ){?>TT 2 no) 2 ^ the Fermi electron temper- 
ature. It turns out that lVr e = ZfTp/Ti ^> 1, since 
in quantum plasmas we usually have Tp > Ti. The dy- 
namics of degenerate electron fluids is governed by the 
continuity equation 



is the electron Fermi speed and tq 
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the momentum equation [2£ 
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the Wigner-Seitz radius, and the sum of the electron 
exchange and electron correlations potential is [HJ V xc — 
0.985(e 2 /e)n 1 / 3 [f + (0.034/a B n 1/3 )ln(I + 18. 37 a s n 1/3 )] , 
where a B = e?i 2 /m*e 2 represents the effective Bohr 
radius. The quantum Bohm potential is [30l - [32j 
Vb = (fi 2 /2m„)(l/y/n)S7 2 y/n. We have thus retained 
the desired quantum forces that act on degenerate 
electrons in our non-relativistic quantum plasma. The 
quantum hydrodynamic equations (l)-(3) are valid 
[26l |34[ |35| if the plasmonic energy density huj pe is 
smaller (or comparable) than (with) the Fermi electron 
kinetic energy k B Tp, where ui pe — (47rnoe 2 /em*) 1 / 2 
is the electron plasma frequency, and the electron-ion 
collision relaxation time is greater than the electron 
plasma period. 

Letting n = no + m, where n\ <«o, we linearize the 
resultant equations to obtain the electron density pertur- 
bation n\ that can be inserted into Eq. (3) . The Fourier 
transformation in space leads to the electric potential 
around an isolated ion. In the linear approximation, we 
have 0, [3 
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where r denotes the position relative to the instantaneous 
position of the point test charge, and the dielectric con- 
stant for a dense quantum plasma with quasi-stationary 
density perturbations is given by 



D = I 



pe 



k 2 (v 2 /3 + v 2 x ) + h'k 4 /4m 2 



(5) 



Here we have denoted v ex = (0.328e 2 /m*er ) 1 / 2 [I + 
0.62/(1 + l8.36a B n /3 )]^ 2 . 

The inverse dielectric constant can be written as 



1_ _ (k 2 /k 2 s ) + qfc 4 /fcf 
D ~ l + itf/kD+ak^/kj 



(6) 



where k s = cj pe / \/v 2 /3 + v 2 x is the inverse 
Thomas-Fermi screening length, and a = 
U 2 ijj 2 e /Arn 2 (v 2 /2> + v 2 x ) 2 measures the importance 
of the quantum recoil effect. If = m e , e = 1 then a 
depends only on r /ag with a B — h 2 /m e e 2 . Inserting 
Eq. ^ into Eq. the latter can be written as 



and Poisson's equation 
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(n - n ) - 47rQ<5(r), 
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where n is the electron number density, u the electron 
fluid velocity, <f> the electric potential, and e the relative 
dielectric permeability of the material (for example, 
for semiconductor quantum wells we have e = 13). 
We have denoted the quantum statistical pressure 
P = (n TO*u 2 /5)(n/n ) 5 / 3 , where u* = hiiir 2 ) 1 / 3 /m»r 



(r) = -2- 

[ 1 4n 2 



(1 + 6) , (1-6) 



k 2 



exp(zk • r)d 3 k, (7) 



where 6 = 1/^1 - 4a, and k 2 ± = k 2 [l + yjl - 4a] /2a. 
Here, we use y/1 — 4a — i\j4a— 1 for a > 1/4. The 
integral in Eq. ([7]) can be evaluated using the general 
formula 



exp(ik ■ r) 
k 2 + k\ 



(Pk = 2tt 



2 exp(— k±r) 



(8) 
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FIG. 1: The value of a as a function of as/ro- The critical 
value a = 1/4 is indicated with a dotted line. 



which involves the Friedel oscillations [38| arising from 
the Kohn anomaly [39| related with the discontinuous 
Fermi surface. 
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FIG. 2: The electric potential as a function of r for a = 
0.625 (dashed curve), a = 1/4 (solid curve) and a = (dotted 
curve). The value 0.627 is the maximum possible value of a 
in our model, obtained for a_g/ro ~ 0.15. 



where the branches of k± must be chosen with positive 
real parts so that the boundary condition <j) — > at r — > 
oo is fulfilled. 

First, for a > 1/4, the solution of the equation 
k 2 ± = k 2 s (lTW4:a-l)/2a yields k± = {k s /Via)[{V^a+ 
1) 1 ^ 2 T«('\/4ck— l) 1 ^ 2 ] = k r ^fiki, and the electric potential 

4>(r) — — [ cos(fcir) + 6* sin(fcir)] exp(— k r r), (9) 

where 6* = l/v^a - 1. In the limit ?i 2 fc 2 > 4m 2 (u 2 /3 + 
w 2 ^ ) , we recover the exponential cosine-screened Coulomb 
potential (37j 

0(r) = — cos(fc s f) exp(— k s r), (10) 
r 

where f = r/(4a) 1 / 4 . Second, for a — > 1/4, we have 
k + = k- = y2k s , and 

#r) = ^(l + ^W(-V2fc.r). (11) 



Third, for cv < 1/4, the expression y/l — 4a is real, and 
we obtain k± = k s (l =p ^1 — 4a) 1 / 2 /v / 2a. The resultant 
electric potential is 

0(r) = ^1" (1 + 6) exp(-fc+r) + (1 - b) exp(-fc_r)l . (12) 
Zr 

We note that in the limit a — > 0, we recover from (fT2|) 
the modified Thomas-Fermi screened Coulomb potential 
4>(r) = (Q/r) exp(— k s r). Furthermore, the newly found 
electrical potential, given by Eq. (12), is significa ntly 
different from that potential [e.g. Eq. (13) in Ref. [36j j . 
indicating that the electric potential is proportional to 
r~ 3 cos(2k F r), where k F — pf/^ — |k||2 is the Fermi 
wave number and p F the Fermi electron momentum] 
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FIG. 3: a) The distance r — d from the test ion charge where 
d<f>/dr = and the electric potential has its minimum, and 
b) the values of the potential cf> and c) its second derivative 
d 2 (f>/dr 2 at r = d. 

Figure 1 shows the value of a as a function of a^/ro, 
with a maximum a at a^/ro ~ 0.15, corresponding 
to a number density no ~ 2 x 10 22 cm -3 (with as — 
5.3 x 10 _9 cm), a few times below solid densities. The 
value of a is above the critical value 1 /4 only for a lim- 
ited range of denities. In Fig. 2, we display the profiles of 
the potential [given by Eq. ©, HU) and ((HJ) for a > 1/4, 
a = 1/4 and a < 1/4, respectively] for different values 
of a. We clearly see the new short-range attractive elec- 
tric potential that resembles the LJ-type potential for 
a > 1/4, while for smaller values of a, the attractive po- 
tential vanishes. Figure 3(a) depicts the distance r = d 
from the test ion charge where d<f>/dr — and the elec- 
tric potential has its minimum, and Figs. 3(b) and 3(c) 
show the values of cj) and d 2 4>/dr 2 at r = d. The value of 
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(d 2 (j)/dr 2 ) determines the oscillation frequency and dis- 
persive properties of the plasma as shown below. For 
a < 0.25, the electric potential and its second deriva- 
tive vanish, and there is no attractive potential associ- 
ated with the stationary test ion charge. Furthermore, 
we note that the shielding of a moving test charge and 
bound states near a moving charge in a quantum plasma 
have been investigated by Else el al. [40j by using the 
Lindhard dielectric function [4l| that ignores the electron 
exchange and electron correlation effects. 

The interaction potential energy between two dressed 
ions with charges Qi and Qj at the position and Vj 
can be represented as Uij(R*ij) = Qj<fii(Rij), where fa 
is the potential around particle i, and = — Vj. For 
a > 1/4, it reads, using Eq.(j9|), 



UijiRij) = ^|exp(-fc r |Ry|) 
x [cos(fcj|Ry|) +6*sin(fe i |Ri i |)] 



(13) 



On account of the interaction potential energy, ions 
would suffer vertical oscillations around their equilibrium 
position. The vertical vibrations of ions in a crystallized 
ion string in quantum plasmas are governed by 



M 



d 2 Szj(t) 
dt 2 



dUij(ri 



(14) 



where dzj(t)[= Zj(t) — Zj ) is the vertical displacement of 
the jth ion from its equilibrium position Zjo, and M the 
mass of the ion. Assuming that 8zj (t) is proportional to 
exp[— i(u;t — jka)], where u and k are the frequency and 
wave number of the ion lattice oscillations, respectively, 
and that Qi = Qj = Q, Eq. (fl4")l for the nearest- neighbor 
ion interactions gives 



(15) 



where S = [2(1 + k r d) + (k 2 - fe?)d 2 ](cos£ + 6* sin£) + 
2fcjd(l + fc r d)(sin£ — 6*cos£). £ = kid, and d is the 
separation between two consecutive ions. We note that 
Eq. (fl~5|) can also be obtained from the formula [HJ 
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M 



d 2 W{r) 



dr 2 



r—d 



(16) 



where the inter-ion potential energy is represented as 
W(r) = (Q 2 /r) exp(— k r r) [cos(fc;r) + 6* sin(fcjr)] for a > 



1/4. Hence, the lattice wave frequency is proportional to 
[d 2 far)/dr 2 }l f J d [cf. Fig. 3(c)], and vanishes for a < 0.25. 

Summing up, we have discovered a new attractive force 
between two ions that are shielded by degenerate elec- 
trons in an unmagnetized quantum plasma. There are 
several consequences of our newly found short-range at- 
tractive force at quantum scales. For example, due to the 
trapping of ions in the negative part of the exponential 
oscillating-screened Coulomb potential, there will arise 
ordered ion structures depending the electron density 
concentration, which in fact controls the Wigner-Seitz 
radius ro. The formation of ion clusters/ion atoms will 
emerge as new features that are attributed to the new 
electric potential we have found here. Finally, under the 
action of the attractive force, we can have the formation 
of Coulombic ion lattices (Coulomb ion crystallization) 
and ion lattice vibrations, as well as the phenomena of 
phase separations at nanoscales in dense quantum plas- 
mas (e.g. from solid to liquid-vapor phases) depending 
upon how one controls the ratio ro/as- Thus, the ratio 
between the interaction energy between the two nearest- 
neighbor ions in the presence of the oscillating exponen- 
tial Coulomb potential and the ion thermal energies, as 
well as the inter-particle spacing are the key parameters 
which will determine a critical point that is required for 
phase transitions in quantum plasmas. In conclusion, the 
present investigation, which has revealed the new physics 
of collective electron interactions at nanoscales, will open 
a new window for research in one of the modern areas of 
physics dealing with strongly coupled degenerate elec- 
trons and non-degenerate mildly coupled ions in dense 
plasmas that share knowledge with cooperative phenom- 
ena (e.g. the formation of ion lattices) in condensed mat- 
ter physics and in astrophysics. Thus, the present inves- 
tigation contributes to enhancing the existing knowledge 
of Wi gne r cr y stallization in two-component Coulomb sys- 
tems [43|, |44j that do not account for an attractive force 
between like charged particles. 
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